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FOREWORD 


Tha  conditions  for  transition  from  primarily  conductive  to  convective  heat 
transfer  in  a  simple  porous  material  have  been  analyzed,  assuming  the  validity 
of  Darcy's  law.  A  stability  criterion,  which  may  apply  in  the  Def lagration-to- 
Detonation  Transition  (DDT)  of  porous  explosives  and  propellants,  has  been 
derived  and  discussed.  This  work  should  be  of  interest  to  those  concerned  with 
propellant  safety  and  DDT  in  porous  reactive  materials. 
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I.  INTRODUCTION 


In  1948  Lapwoodl  analyzed  Che  heat  flow  in  an  infinite,  horizontal, 
fluid-saturated,  homogeneous,  porous  slab.  By  adopting  the  mathematical 
technique  used  by  Lord  Rayleigh  in  1916^  in  an  analogous  problem  (the 
stability  of  a  liquid  layer  heated  from  below),  Lapwood  showed  that  the  flow  is 
not  always  stable,  but  that  there  exists  the  same  stability  condition  that 
Rayleigh  had  found  in  the  case  of  the  liquid  layer.  The  stability  condition 
states  that,  if  the  temperature  gradient  exceeds  a  certain  critical  value,  one 
has  the  onset  of  instabilities  associated  with  fluid  convection. 


This  analysis  was  conducted  assuming  zero  viscosity,  which  drastically 
simplifies  the  problem.  Only  recently,  Kassoy  and  Zebib^>^>5  have 
investigated  the  same  problem  by  taking  into  account  the  viscosity  and  using  the 
more  realistic  configuration  of  a  confine  ’  medium. 


The  onset  of  instabilities  corresponds  to  a  transition  from  a  conductive  to 
convective  heat  transfer  which  manifests  itself  in  a  dramatic  change  in  the 
thermal  conductivity  characteristics  of  the  porous  material.  As  early  as  1933, 
Kannuluik^  showed  that  the  thermal  conductivity  of  powders  shows  a  sharp 
increase  when  they  are  permeated  by  a  gas.  For  example,  for  evacuated  powders 
such  as  carborundum,  the  thermal  conductivity  is  typically  of  the  order  of  C. 1  x 
10“ 5  cal/cm-sec-°C.  However,  when  carborundum  powder  is  filled  with 
hydrogen,  whose  conductivity  is  38.0  x  10“^  jn  the  same  units,  the 
conductivity  of  the  powder  increases  to  300  x  10” 5  cal/cm-sec-°C,  3  x  10^ 
times  the  thermal  conductivity  of  the  evacuated  powder.' 


^Lapwood,  E.  R.,  "Convection  of  a  Fluid  in  a  Porous  Medium,"  Proc.  Carob.  Phil. 
Soc .,  44,  5C8-521,  1948. 

2bord  Rayleigh,  "On  Convection  Currents  in  a  Horizontal  Layer  of  Fluid,  When 
the  Higher  Temperature  is  on  the  Under  Side,"  Phil.  Mag., 6,  32,  529-546,  1916. 

^Kassoy,  D.R.  and  Zebib,  A.,  "Variable  Viscosity  Effects  on  the  Onset  of 
Convection  in  Porous  Media,"  Phvs.  of  Fluids.  12 ,  18,  1649-1651,  1975. 

^Zebib,  A.  and  Kassoy,  D.  R.,  "Onset  of  Natural  Convection  in  a  Box  of  Water- 
Saturated  Porous  Media  with  Large  Temperature  Variation,"  Phys.  of  Fluids,  2^0 
(1),  4-9,  1977 

^Zebib,  A.,  "Onset  of  Natural  Convection  in  a  Cylinder  of  Water-Saturated 
Porous  Media,"  Phys.  Fluids,  21_  (4),  699-700  , 137S. 

6Ka  nnuluik,  W.  G.,  "Conduction  of  Heat  in  Powders,"  Proc.  Poy.  Soc.tA141, 
144-158,  1933.  “ 
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This  phenomenon  can  be  explained  only  by  assuming  that  a  convective  process 
is  established  and  that  the  heat  is  transferred  to  the  grains  of  powder  by 
collision  with  the  fluid  molecules. 

The  effect  of  the  gas  imbedded  in  a  porous  medium  is  illustrated  in  Figure 

I,  which  also  shows  the  importance  of  grain  sizes  and  density  for  carborundum 
powder  (See  Table  I  for  density  and  particle  size  parameters  for  carborundum). 
The  effect  of  different  gases  imbedded  in  the  powder  (hydrogen  and  air)  is  also 
shown.  Figure  1  shows  also  that  the  conductivities  increase  with  pressure, 
although  they  tend  to  reach  a  plateau  at  pressures  of  the  order  of  one 
atmosphere,  so  that  for  high  pressures  one  can  assume  that  the  conductivities 
become  pressure-independent.  This  seems  reasonable  if  one  considers  the  data 
shown  on  curves  of  conductivity  vs  pressure  in  Figure  2  for  two  typical 

gases — namely,  air  and  hydrogen. 

These  results  appear  to  play  an  important  role  in  the  def lagration-to- 
detonation  transition  in  high  explosives  because,  as  is  well  known,  this 
transition  is  preceded  by  another  transition,  that  from  conductive  to  convective 
burning.  An  analysis  to  determine  the  conditions  necessary  for  this  transition 
will  now  be  considered. 

II.  STABILITY  ANALYSIS 

The  xy-plane  contains  a  horizontal  slab  of  a  fluid-saturated  homogenous, 
porous  material  with  thickness  h,  so  that  0  <_  z  h. 

In  absence  of  viscosity,  the  Navier-Stokes  equation  reduces  to  the  Euler 
equation : 

P^  =  -Vp-F  (1) 

Dt 


with 

*  i_  +  (T.  T  (2) 

Dt  3 1 


p ,  the  density,  V,  the  macroscopic  velocity  of  the  element  of  fluid  at  the_^ 
point  (x,y,z),  whose  components  are  u,v,  and  w,  with  p,  the  pressure,  and  F,  the 
system  of  forces  acting^n  the  element  of  fluid.  This  force  consists  of  two 
components,  the  weight  W  =  (0,0,pg),  with  g,  the  acceleration  of  gravity,  and 
a  resistance  component  for  the  porous  material  given  by  Darcy's  law,  so  that  one 
has  for  the  total  force: 


(3) 
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K  being  the  permeability  and  lT,  the  unit  vector  in  the  direction  of  the  z-axis 
The  gravitational  component  pgk,  although  negligible  in  comparison  to  the 
Darcy  term  (which,  in  the  case  of  low  permeability  materials,  such  as 
propellants,  is  large)  has  been  included  for  the  sake  of  completeness. 

To  the  fluidynamic  equation,  one  must  add  the  heat  conduction  equation: 


51  =  k?29 
Dt 


(4 


where  0  is  the  temperature  and  <  the  thermal  diffusivity.  To  Equations  1-4, 
one  must  also  add  two  other  equations,  a  continuity  equation: 


I  2s.  -  -IT .  T  (5 

p  Dt 


and  an  equation  for  thermal  expansion: 


p  =  P  c  (  l-o  0  ) 


(6 


where  pc  is  the  value  of  p  at  0  =0,  and  a  is  the  bulk  thermal 
expansion  coefficient. 


Equations  (l)-(6)  can  be  linearized  by  assuming  that  velocity  and  the 
departures  of  temperature,  density,  and  pressure  from  the  equilibrium  state 
small  so  that  second  order  terms  can  be  assumed  negligible. 


ari 


If  the  variations  in  density  as  well  as  the  (V  .  V)V  term  are  considered 
negligible,  one  has 


p  3jT  =  -  Tp  -  ££  V  -  p  gif  , 
at  k  ' 


(7 


51  =  k?29  , 

(8 

Dt 

r .  t  =  o  , 

(9 

0  -  Or-  (1-Ct  0) 

(10 
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The  subscript  "0"  labels  the  parameters  in  the  equilibrium  state.  The 
assumed  perturbed  from  the  equilibrium  state  and  the  perturbations  are 
usual,  considered  small.  We  have: 


system  is 
as 


9=eo+0=-fi2+e  , 

(ID 

p  =  Po  +  P  55  ~  P  08z  +  P  » 

(12) 

°o  =  pc  (l“a®0)  , 

(13) 

P  =  pc  1-a  (6  +0  Q)  =  pQ  -  pcae  , 

(14) 

where  (1  is  the  temperature  gradient. 

Since  the  term  V  .  V0  is  negligible,  equation  (8)  reads 


-  8w  +  =  kV  0  . 

3  t 


(15) 


By  taking  the  divergence  of  (7),  since  (10)  holds,  one  has 


2 

V  P 


PcStt 


ae_ 

3  z 


0. 


(16) 


The  third  component  of  (7)  gives 


-^»-°os*epc«o  , 


(17) 


and  with  (15)  and  its  time  derivative,  w  can  be  eliminated,  and  one  h 


as 


-  l£+  p  cga0  = 


3  z 


i  °'(tt  *  t)  (It  -  "2) 9  *  'os 


(18) 


By  taking  the  Laplacian  of  (18)  and  the  z  derivative  of  (16),  the  pressure  also 
can  be  eliminated,  and  one  readily  has  the  fourth  order  partial  differential 
equation  for  the  temperature  0 , 
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2  2  2 
7  8  *  ga  (5  7  0  , 


where  7^  ,  is  the  partial  Laplacian,  namely 


2  2  2 
7  =  3 _  +  _3 _ 

1  3  x2  3  y2 


Let  us  now  consider  a  simple  case,  with  the  following  boundary  conditions: 


0=0  at  z  =  0  and  z  =  h 


JLJL  =0  at  z  =  0  and  z  =  h 


The  search  for  a  solution  of  Equation  (19)  of  the  form 


0  =  R(x,y,t)  Z(z) 


imposes  on  z,  through  the  boundary  conditions  (21)  and  (22),  the  following 
conditions , 


Z(0)  =  Z(h)  =  0  , 


Z"(0)  =  Z"(h)  =  0  , 


so  that  by  taking 


Z ( z )  =  A  sin  uz 


one  has 


sinwh  =  0  i.e.,  uh  *  nr 


Now  we  try  a  solution  of  the  form 
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0 


AsinX xsiny ysin  I2L 
h 


vt 

z  e 


(28) 


By  substitution  into  Equation  (19),  one  has 


(  2  2  2  1 

/  2 

2  2’ 

\  2 

♦i) 

)v  +  £_  (a  +  n  it  )( 

j  *2  j 

+  n  it 

j  =  gafia 

where 

a^  =  (y  2  +  x  2)  > 

a  =  £  , 

1  K 


(29) 


C 

(2 


a  2 


*  (a  +  n  ir  )  , 
h2 


(32) 


2 

a  =  £”3  (j  (33) 

3  a2  +  n42 


and  Equation  29  reads 


(v  +  aj)  (v  +  02)  =  013 


(34) 


which  admits  the  following  roots 


v 


1/2  (ox  +  o2)  +  1/2 


J  2 

^f(a^-a2)  +  403 


(35) 


When  v  is  complex,  its  real  part,  -1/2(03  +  a2),  is  always  negative 
and  the  perturbation  decays  with  time,  since  03  and  a2  are  positive.  If 
we  require  the  perturbation  to  persist  then  the  roots  must  be  real  and  positive, 
i.e. , 
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2 

(aj  -  02)  +  4ct3>0, 


(36) 


and 


+  <*2)  + 


- 1 - 

02)  +  4a  3  >  0  • 


(37) 


If  (3  >  0,  (35)  is  certainly  satisfied,  and  the  roots  are  real. 

The  value  of  8  s  13cr  which  results  in  a  value  of  zero  for  v  is  the  one 
for  which,  when  13  <  I3cr,  the  perturbation  dies  down,  and  when  (3  >  I3cr> 
it  increases  exponentially.  This  value  is  found  by  solving  the  equation: 


2  2 
(01-02)  +403  =  (01+02)  , 

(a  +  a  )2  -  (a  -  a,)2 

a  =  I _ t _ i _ , 

3  4 

« 3  =  ala2  . 


which  gives 


,  2  2  2,2 

(a  +  n  tt  ) 


c  r 


Koh' 


(38) 


(39) 


(40) 


(41) 


The  minimal  value  of  (3cr  =  (3cr  is  obtained,  since  n  is  an  integer, 

2  2 

for  n  =  1  and  for  a  value  of  a,  such  that  the  function  3  +  w  is  minimum. 

a 

This  minimum  is  readily  recognized  for  a2  =  tt  2 ,  so  that 


2 

B  "  *L4  ,  («) 

cr  Koh2 


and  for 
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2 

(5  <  4tt  k 
Kah2 


(43) 


the  flow  is  stable  and  the  perturbation  rapidly  dies  exponentially.  For 


2 

13  >  4tt  k  f 
Kah2 


(44) 


the  perturbation  increases  exponentially  and  the  flow  is  unstable.  This 
represents  the  condition  for  the  transition  from  conductive  to  convective  heat 
transfer  in  a  fluid-filled  porous  material  in  the  hypothesis  that  Darcy's  law 
holds . 

I  IT.  THE  ROLE  OF  THE  STABILITY  CONDITION  IN  THE  BURNING  OF  POROUS  MATERIALS  AND 
PROPELLANTS 


The  stability  condition  discussed  in  the  previous  pages  appears  to  be  of 
particular  interest  in  the  study  of  the  combustion  of  porous  material  and 
propellants.  It  involves  all  the  typical  parameters  one  expects  to  encounter  in 
the  process;  namely,  the  temperature  gradient,  the  thermal  diffusivity,  the 
thermal  expansion  coefficient  of  the  fluid  filling  the  pores,  the  permeability 
of  the  porous  material,  and  the  thickness  of  the  material  through  which  the 
thermal  gradient  is  established. 

The  stability  condition  states  that  when  the  temperature  gradient  (3  is  such 

that 


2 

13  >  *  , 

Kah2 


the  heat  transfer  dramatically  changes  character  from  conductive  to  primarily 
convective.  Analogous  criteria  are  expected  to  apply  to  the  burning  of  modern 
solid  propellants  with  a  high  content  of  solid  materials.  In  this  case,  during 
combustion,  ahead  of  the  flame,  there  exists  a  layer  of  propellant-heated, 
although  not  yet  burning.  Although  the  temperature  is  initially  insufficient  to 
ignite  the  propellant,  it  increases  with  time  and  can  reach  a  value  at  which  the 
solvent  and  the  binder  begin  to  evaporate  or  decompose  so  that  a  porous  layer 
whose  thickness  varies  with  time,  begins  to  form  ahead  of  the  flame.  The 
previous  considerations  then  apply.  In  the  case  of  small  propellant  grains,  the 
total  burning  time  is  probably  too  short  for  such  a  phenomenon  to  occur.  The 
advent  of  longer  range  solid-propellant  requires  longer  burning  propellant 
grain,  so  that  the  process  described  above  could  develop. 
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IV.  AN  ESTIMATE  OF  THE  THERMAL  GRADIENT  ESTABLISHED  IN  A  PROPELLANT  GRAIN 


The  problem  of  heat  transfer  in  thermal  conductors  when  the  heat  sources  are 
moving,  is  a  problem  of  foremost  importance  in  several  fields,  such  as  metal 
welding,  and  many  studies  have  been  made  on  the  subject.2*®  However,  since, 
from  a  mathematical  point  of  view,  the  rigorous  treatment  is  of  great 
difficulty,  this  analysis  has  been  confined  to  the  quasi-stationary  regime. 

This  regime  is  experimentally  verified^  and  is  seen  in  the  fact  that,  if  the 
sample  is  long  enough,  the  temperature  ahead  of  the  moving  source  rapidly 
approaches  a  constant  value.  In  different  terms,  an  observer  who  is  moving  with 
the  source  does  not  detect,  after  a  certain  time,  any  change  in  the  temperature 
distribution  in  the  sample  ahead  of  the  source.  While  the  quasi-stationary 
solution  is  often  of  practical  utility,  the  transient  regime  is  important  in 
many  cases. 

For  the  sake  of  simplicity,  we  will  make  an  assumption,  which  will  be 
removed  in  a  successive  study,  that  the  sample  is  non-reactive.  The  burning 
process  is  then  schematized  by  substituting  for  the  flame  front  a  heat  source 
which  moves  along  the  x-axis  with  a  constant  speed  v  and  a  constant  heat  rate 
q.  The  evaluation  of  the  temperature  distribution  T(x,y,z,t)  will  be  determined 
by  the  heat  conduction  equation: 


3  t  -f  3  T  +  3  T  x  2X  3  T  j 
3x2  3  y2  3  12  3t 


(45) 


where 


(46) 


c  being  the  specific  heat  of  the  solid  material,  o  its  density  and  k  the 
thermal  conductivity  coefficient  of  the  material.  The  transformation 


?  -  x  -  vt, 


(47) 


brings  Equation  (45)  into  the  following  form: 


^Campolattaro,  A.  A.,  Numerical  Modelling  of  the  Burning  of  Highly  Exothermic 
Mixtures,  NSWC  report  1977  (To  be  published). 

^Rosenthal,  D.,  The  Theory  of  Moving  Sources  of  Heat  and  Its  Applications  to 
Metal  Treatment,  Trans.  Amer,  Soc.  Mech  Engrs.,68  849-866,  3.946,  and  included 
bibliography. 

^Bornefeld,  H. ,  Temperature  Measurement  in  Fusion  Welding,  Technische 
Zentralblatt  Fur  Praktische  Metal  Bearbeitung  Vol.  43_.  pp.  14-18,  1933. 
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V  T  +  2A  v  ±1 
3C 


2A  11  =  0. 
3  t 


(48) 


where  the  Laplacian  is  the  Eulerian  one,  i.e.  the  operator  Laplacian  involving 
the  coordinates  system  moving  with  the  heat  source,  i.e. 


2  ,2  ,2  o2 

V  =i _  +i _  +  1 _  .  (49) 

3? 2  3y2  3  z2 


By  putting 


T  =  T0  +  e~<*Z  (C  ,y,z,t)  (50) 


where  T0  is  the  initial  temperature,  Equation  (48)  becomes 


2 

7  <t  +  2Uv-a )  li 
3C 


2A  1±  +  a  (a -2A  v)  4>=0. 
3  t 


(51) 


and  by  taking 


a  =  A  v , 


(52) 


the  term  11  can  be  eliminated  and  one  has 
3C 


2 

V  4> 


2A  ii  -  (Av) 
3  t 


$=0. 


(53) 


3y  putting 


2 

-A  v  t 

$  =  e  2  vt»  , 


(54) 


Equation  (53)  reduces  further  to  the  ordinary  heat  conduction  equation,  i.e. 
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r 


s 

1 


2 

7  \p 


-  2A  j±  =  0  , 
3t 


(55) 


The  solution  of  (48)  which  satisfies  the  boundary  condition  that  at  the  flame 
(C  =0),  the  temperature  remains  constant  and  equal  to  Tf,  for  the 
one-dimensional  case,  will  be 


T  -  T0 


-Av£  2 

e  e  ip(£,t)  , 


(56) 


so  that  the  initial  condition  for  ip  is 


<P  (C  ,  t)  =  0  for  t  =  0, 


and  at  the  boundary, 


(57) 


-A  v‘ 


Tf  -  T0  =  e  'P  (0,  t )  , 


(58) 


Equation  (55)  must  be  integrated  with  the  initial  and  boundary  conditions  given 
by 


<P  (£  ,  t)  =  0  for  t  =  0, 


Av 


(59) 


<p(£,t)  =  (Tf  -  T0)  e  for  £  =  0. 


This  integral  is  readily  found^ 


(60) 


l^Carslaw,  H.  S .,  and  Jaeger,  J.  C.,  Conduction  of  Heat  in  Solids,  Oxford  at  the 
Clarendon  Press,  II  Ed., 1959,  page  63. 
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XvC  2 

e  e  (T  -  Tc)  - 


2  2 
X  v  t '  -  X  £  y 

2  2  C  t-t  *  7 

e _ e _ d_t 

(t  -  t')3/2 


(61) 


By  putting 

t  -  t '  =  \i  ,  ( 62) 

Equation  (61)  reads 


T-T 


o 


jn  (T  -  T  )e  e 
1  2ir  f  O 


/w-3/ 2 

oJ 


2  2 
-Xvu  -  X£ 

e  2  e  2y  ay 


(63) 


or 


T 


T 

o 


(T 

f 


T  )C 
o 


-X  vC 


2  2  2 

t  -  X  C  +  V  ii 


(64) 


which  gives  the  temperature  at  any  time  t  and  at  a  distance  5  from  the 
propagating  flame. 

An  easy  check  of  this  solution  can  be  made  by  set  ’ing  v“0  in  Equation  (64) 
to  determine  whether  or  not  it  reduces  to  the  solution  for  a  semi-infinite  slab 
with  the  boundary  at  £=0  kept  at  constant  temperature  Tf  and  the  initial 
temperature  constant  and  equal  to  T0,  which  is  a  classical  problem. 

For  v=0,  Equation  (64)  reads: 


^ ‘Cars law,  H, 


C. ,  Conduction  of  Heat  in  Solids. 
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-  i  i_ 

f  -3/2  2  V 

re)c  J  y  e  dy  , 

o' 


which,  by  setting 


readily  gives: 


T  -  T„  - 


=  (Tf  -  T0)  erfc  C  JII 


which  is  indeed  the  solution  given  by  Carslaw  and  Jaeger. 

It  is  also  useful  to  evaluate  the  steady-state  solution  of  Equation  (48). 
Neglecting  the  time  derivative  in  Equation  (48),  we  have  as  the  steady-state 
equation, 


7  T  +  2Av  II  =  0. 


By  putting 


T  =  T0  +  e~x^  # 


Equation  (65)  becomes 


7  2$  -  (Xv)2  4>  =  0, 


which  reduces  to,  in  the  one-dimensional  case: 


-  (A v)2  $  =  0, 

dC2 


Carslaw,  H.  S.  and  Jaeger,  J.  C. ,  Conduction  of  Heat  in  Solids,  p.  60. 
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whose  general  integral  is 


$  =  AeXv?  +  Be“Xv?  ,  (72) 

with  A  and  B  as  arbitrary  constants  of  integration. 

The  boundary  conditions  for  T  which  are  of  physical  interest  are: 


T  =  Tf  for  C  =  0,  and  lim  T  =  T0, 

£-►00 


(73) 


and,  from  (69)  and  (72),  one  readily  has 


T  =  Tq  +  (Tf  -  Tq)  e-2vA£. 


(74) 


If  Tr  is  the  temperature  at  which  some  chemical  reaction  or  phase  change 
occurrs  within  one  or  more  of  the  components  of  a  mixture  undergoing  combustion, 
then  the  distance  Cr,  ahead  of  the  flame,  at  which  the  burning  material  has 
reached  temperature  Tr  is  readily  evaluated  from  (74),  and  one  has 


This  distance  is  only  of  interest  if  we  have 


Cr  >  vr 


(76) 


where  x  is  the  reaction  time  at  the  temperature  Tr.  Equation  (76)  combined 
with  Equation  (75)  gives  the  following  relation  among  the  physical  parameters  of 
the  system  for  a  reaction  or  phase  change  to  occur: 


T  T 

A v^x  <  In  £  0 

Tr  -  T0 


(77) 


Equation  (64)  may  be  integrated  numerically  to  evaluate  the  temperature  at 
any  time  x  at  a  distance  C  from  the  propagating  flame  for  materials  whose 
burn  velocity  (v),  flame  temperature  (Tf)  and  thermal  diffusivity  are  known. 
These  parameters  are  summarized  for  three  materials  in  Table  2.  Equation  (64) 
was  integrated  for  each  material  at  various  distances  (£ )  ahead  of  the  flame 
as  a  function  of  time.  The  results  are  shown  in  Figures  3(a),  3(b)  and  3(c). 

The  time  required  to  establish  a  "steady-state"  temperature  varies  directly  with 
the  distance  ahead  of  the  flame  (f )  and  the  thermal  diffusivity,  as  expecced. 
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Calculations  were  also  made  to  determine  "steady-state”  temperature  rise  as  a 
function  of  distance  ahead  of  the  flame  front  by  integrating  at  each  of  10 
distances  numerically  until  convergence  is  reached.  The  results  of  these 
calculations  for  materials  listed  in  Table  1  are  shown  in  Figures  4(a),  4(b), 
and  4(c).  The  calculated  temperature  as  C^O,  though  not  shown  in  these 
figures,  does  indeed  approach  the  flame  temperature.  From  these  calculations, 
an  estimate  can  be  made  of  the  thermal  gradient  existing  in  porous  materials 
burning  conductively  to  assess  when  the  transition  conditions  for  the  onset  of 
convective  flow  are  reached.  Measurements  of  the  permeability  (K)  of  these 
propellants  as  well  as  the  thermal  expansion  coefficient  of  their  combustion 
products  (a)  must  be  made  in  order  to  check  the  validity  of  the  model. 

V.  CONCLUSIONS 

A  condition  has  been  found  which  governs  the  transition  from  conductive  to 
convective  burning  in  nonreactive  analogs  of  porous  propellants  and  explosives. 
The  relevant  physical  parameters  in  the  transition  are,  as  one  would  expect,  the 
temperature  gradient  in  the  material,  its  thermal  diffusivity,  and  the  thermal 
expansion  coefficient  of  its  combustion  gases.  It  hj?s  been  previously 
well-established1^’  4  that  the  conductive-convective  flow  transition  precedes 
that  from  deflagration-to-detonation.  A  mechanism  has  been  proposed  whereby  an 
initially  non-porous  propellant  grain  may  become  porous  through  conductive 
heating  ahead  of  the  flame  and  subsequent  phase  change  or  thermal  degradation  of 
propellant  ingredients.  The  establishment  of  a  sufficiently  high  thermal 
gradient  throughout  this  newly-formed  porous  layer,  whose  thickness  increases 
with  time,  then  defines  the  condition  for  the  onset  of  primarily-convective  flow. 

Experimental  work  will  be  done  in  the  near  future  to  obtain  the  physical 
parameters  necessary  to  verify  the  validity  of  the  condition  for 
conductive-convective  flow  transition  in  porous  materials  and  its  role  as  a 
precursor  to  DDT. 
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Griffith,  N.,  and  Gloocock,  J.  M.,  "The  Burning  to  Detonation  of  Solid 
Explosives,"  J.  Chem.  Soc.,4154,  I960  . 

^Bernecker,  R.  R.,  and  Price,  D.,  "Studies  in  the  Transition  from 

Deflagration  To  Detonation  in  Granular  Explosives-I  Experimental  Arrangement 
and  Behavior  of  Explosives  which  Fail  To  Exhibit  Detonation,"  Combust.  Flame, 
22,  111-117,  1974. 
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FIGURE  1.  THERMAL  CONDUCTIVITY  (k)  CURVES  FOR  CARBORUNDUM 
WITH  H2,  SHOWN  IN  DOTTED  LINE,  AND  FILLED  WITH  AIR, 
SHOWN  IN  SOLID  LINE,  VS  PRESSURE  (REF.  6) 
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FIGURE  2.  CONDUCTIVITY  OF  HYDROGEN  AND  AIR  AS  A  FUNCTION 
OF  PRESSURE  (REF.  6) 


ABOVE 


NSWC 


FIGURE  3c.  CALCULATED  TEMPERATURE 
DISTANCES  AHEAD  OF  THE  FI 
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DISTANCE  AHEAD  OF  FLAME  FRONT  (CM) 


FIGURE  4a.  CALCULATED  “STEADY-STATE"  TEMPERATURE  RISE  AS 
A  FUNCTION  OF  DISTANCE  AHEAD  OF  THE  FLAME  FRONT 
FOR  EJC  PROPELLANT  (AT  t  =  0.045  sec) 


TEMPERATURE  ABOVE  AMBIENT  (DEGREES  KELVIN) 
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0  0.004  0.008  0.012  0.016  0.20 

DISTANCE  AHEAD  OF  FLAME  FRONT  (CM) 


FIGURE  4c.  CALCULATED  "STEADY-  STATE"  TEMPERATURE  RISE  AS 
A  FUNCTION  OF  DISTANCE  AHEAD  OF  THE  FLAME  FRONT 
FOR  PYRONOL  PYROTECHNIC  (AT  t  =  0.006  sec) 
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TABLE  1.  DENSITY  OF  SOLID  MATERIALS  (p)  AS  WELL  AS  POWDERED 
MATERIALS  <p0)  FOR  DIFFERENT  TYPES  OF  CARBORUNDUM. 

THE  AVERAGE  LINEAR  DIMENSION  (I)  OF  THE  POWDERED  MATERIAL 
IS  ALSO  SHOWN.  (REF.  6) 


CARBORUNDUM 

IjjfipPH 

no.  go 

NO.  280 

NO.  600 

p,  g  cm'3 

1  80 

1  •  89 

1  -84 

1  •  54 

Pa.  g  cm'3 

3  -  20 

3-20 

3  •  20 

3  20 

1,  cm 

0  055 

0  •  0194 

0  -  0061 

0  •  0027 

TABLE  2.  SELECTED  PROPERTIES  OF  REACTIVE  MATERIALS  USED 
FOR  CALCULATION 


MATERIAL 

BURN 

VELOCITY 

(cm/sec) 

FLAME 

TEMPERATURE 

THERMAL 

DIFFUSIVITY 

(cm2/sec) 

EJC  PROPELLANT 

1.372 

3913.5 

375.3 

TP-H-1148  PROPELLANT 

1.067 

3429.6 

229.6 

PYRONOL  PYROTECHNIC 

41.58 

2473.0 

2.58 
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